Can all cosmological observations be accurately interpreted with a unique geometry? 
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The recent analysis of the Planck resuhs reveals a tension between the best fits for {flmOjHo) 
derived from the cosmic microwave background or baryonic acoustic oscillations on the one hand, 
and the Hubble diagram on the other hand. These observations probe the universe on very different 
scales since they involve light beams of very different angular sizes, hence the tension between 
them may indicate that they should not be interpreted the same way. More precisely, this letter 
questions the accuracy of using only the (perturbed) Friedmann-Lemaitre geometry to interpret all 
the cosmological observations, regardless of their angular or spatial resolution. We show that using 
an inhomogeneous "Swiss-cheese" model to interpret the Hubble diagram allows to reconcile it with 
the Planck results. Such an approach does not require to invoke new physics nor to violate the 
Copernican principle. 



PACS numbers: 98.80.-k, 04.20.-q, 42.15.-i. 



The standard interpretation of cosmological data re- 
lies on the description of the universe by a spatially ho- 
mogeneous and isotropic spacetime with a Friedmann- 
Lemaitre (FL) geometry, allowing for perturbations [T]. 
The emergence of a dark sector, including dark matter 
and dark energy, emphasizes the need for extra-degrees 
of freedom, either physical (new fundamental fields or in- 
teractions) or geometrical (e.g. a cosmological solution 
with lower symmetry). This has driven a large activity 
to test the hypotheses |2j of the cosmological model, such 
as general relativity or the Copernican principle. 

The recent Planck data were analyzed in such a frame- 
work P] in which the cosmic microwave background 
(CMB) anisotropics are treated as perturbations around 
a FL universe, with most of the analysis performed at lin- 
ear order. Non-linear effects remain small [4] and below 
the constraints on non-Gaussianity derived by Planck [5]. 
The results nicely confirm the standard cosmological 
model of a spatially Euclidean FL universe with a cos- 
mological constant, dark matter and initial perturbations 
compatible with the predictions of inflation. 

Among the constraints derived from the CMB, the 
Hubble parameter Hq and the matter density parame- 
ter J7mO are mostly constrained through the combina- 
tion fimo'i^j where Hq ^ lOOh km/s/Mpc. It is set by 
the acoustic scale 6** = Vs/Da, defined as the ratio be- 
tween the sound horizon and the angular distance at the 
time of last-scattering. The measurement of seven acous- 
tic peaks allows to determine 9^, with a precision better 
than 0.1%. The constraints on the plane (ilmOi^o) sue 
presented in Fig. 3 of Ref. [3^ and clearly show this de- 
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generacy. The marginalized constraints on the two pa- 
rameters were then derived 3\ to be 

Ho = (67.3±1.2) km/s/Mpc, r^^o = 0.315±0.017 (1) 

at a 68% confidence level. It was pointed out (see 
§ 5.2-5.4 of Ref. [3]) that the values of Hq and 17^0 
are respectively low and high compared with their val- 
ues inferred from the Hubble diagram. Such a trend 
was already indicated by WMAP-9 [6] which concluded 
Ho = (70 ± 2.2) km/s/Mpc. 

Regarding the Hubble constant, two astrophysical 
measurements are in remarkable agreement. First, the 
estimation based on the distance ladder calibrated by 
three different techniques (masers. Milky way cepheids 
and Large Magelanic cloud cepheids) gives ^ Hq = 
(74.3 ± 1.5 ± 2.1) km/s/Mpc, respectively with statis- 
tical and systematic errors. This improves the ear- 
lier constraint obtained by the HST Key program [5], 
Ho = (72 ± 8) km/s/Mpc. Second, the Hubble dia- 
gram of type la supernovae (SNe la) calibrated with 
the HST observations of cepheids leads [S] to Hq = 
(73.8±2.4) km/s/Mpc. Other determinations of the Hub- 
ble constant, e.g. from VLBI observations [10] or from 
the combination of Sunyaev-Zel'dovich effect and X-ray 
observations [TT] , have larger error bars and are compat- 
ible with both the CMB and distance measurements. 

Besides, the analysis of the Hubble diagram of SNe la 
leads to a lower value of f^mo - e.g. 0.222±0.034 with the 
SNLS 3 dataset [12] - compared to the constraint ^ by 
Planck. As concluded in Ref. [2| , there is no direct incon- 
sistency, and it was pointed out that there could be resid- 
ual systematics not properly accounted for in the SNe 
data. Still, it was stated that "the tension between CMB- 
based estimates and the astrophysical measurements of 
Ho is intriguing and merits further discussion" . 

Interestingly, the CMB constraints on (Omo,^o) are 
in excellent agreement with baryon acoustic oscillation 



(BAO) measurements [T3], that allow to determine the 
angular distance up to redshifts of order 0.7. The com- 
mon point between the CMB and BAO measurements 
is that they involve light beams much larger than those 
involved in astronomical observations. Indeed, a pixel 
of Planck'' s high-resolution CMB maps corresponds to 
5 arcmin [2], while the typical angular size of a SN is 
10~^ arcsec. This means that the two kinds of observa- 
tions probe the universe at very different scales. More- 
over, for both the CMB and BAO measurements the cru- 
cial information is encoded in correlations, while SNe ob- 
servations rely on "1-point measurements" (we are inter- 
ested in the luminosity/redshift of each SN, not in the 
correlations between several SNe). Because of such dis- 
tinctions one can expect from the two classes of cosmolog- 
ical observations to be affected differently by the inhomo- 
geneity of the universe, through gravitationnal lensing. 

The effect of lensing on CMB measurements is essen- 
tially due to the large-scale structure, and it can be taken 
into account in the framework of cosmological perturba- 
tion theory at linear order [TF (see however Ref. fT6^ for a 
discussion about the impact of strong inhomogeneities) . 
We refer to Ref. [T7] for a description of the lensing ef- 
fects on BAO measurements. Regarding the Hubble di- 
agram, the influence of lensing has also been widely in- 
vestigated [18j . The propagation of light in an inhomo- 
geneous universe gives rise to both distortion and mag- 
nification. In principle, most images are demagnified be- 
cause their lines of sight probe underdense regions, while 
some are amplified due to strong lensing. It shall thus 
induce a dispersion of the luminosities of the sources, 
that is, an extra-scatter in the Hubble diagram [12]. Its 
amplitude can be determined from the perturbation the- 
ory [2D] and substracted [H]. However, a considerable 
fraction of the lensing effects arises from subarcminute 
scales, which are not probed by shear maps smoothed on 
arcminute scales [22) . 

The tension on (JlmOi ^o) rnay indicate that, given the 
accuracy of the observations achieved today, the use of a 
(perturbed) FL geometry to interpret the astrophysical 
data is no longer adapted. More precisely, the question 
that we want to raise is whether the use of a unique space- 
time geometry is relevant for interpreting all the cosmo- 
logical observations, regardless of their angular or spatial 
resolution and of their location (redshift). Indeed, each 
observation is expected to probe the universe smoothed 
on a typical scale related to its resolution, and this can 
lead to fundamentally different geometrical situations. In 
a universe with a discrete distribution of matter, the Rie- 
mann curvature experienced by a beam of test particles 
or photons is dominated by the Weyl tensor. Conversely, 
in a (statistically spatially isotropic) universe smoothed 
on large scales it is dominated by the Ricci tensor. Both 
situations correspond to distinct optical properties [13]. 

In the framework of geometric optics, a light beam is 
described by a bundle of null geodesies. All the infor- 
mation about the size and the shape of a beam can be 
encoded in a 2 x 2 matrix V^^ called the Jacobi map 



(see Ref. [23] for futher details). In particular, the angu- 
lar and luminosity distances read respectively 
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where z denotes the redshift. The evolution of the Ja- 
cobi map with light propagation is governed by the Sachs 
equation [25] [26] 
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where v is an affine parameter along the geodesic bundle. 
The term TZab which controls the evolution of V^^ is a 
projection of the Riemann tensor called the optical tidal 
matrix. It is defined by TZab = Rfivapk'^ k"" s^s^^ where 
k^^ is the wave vector of an arbitrary ray, and the Sachs 
basis {s^} A=i,2 spans a screen on which the observer 
projects the light beam. Because the Riemann tensor 
can be split into a Ricci part and a Weyl part, the optical 
tidal matrix can also be decomposed as 
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Weyl lensing 



with $00 = ~{l/2)R^,kt'k'' and *o = -(l/2)C7^,„^(s^- 
\S2)k'^{si — isj)^^- It clearly appears in Eq. Q that 
the Ricci term tends to isotropically focus the light beam, 
while the Weyl term tends to shear and rotate it. The 
behaviour of a light beam is thus different whether it 
experiences Ricci-dominated lensing (large beams, e.g. 
CMB measurements) or Weyl-dominated lensing (nar- 
row beams, e.g. SNe observations). It has been demon- 
strated [17] that for the lensing caused by a distribution 
of pointlike particles, the Ricci-lensing behaviour is in- 
deed recovered when the size of the beam increases, as 
the light beam encloses more and more particles. 

This Ricci- Weyl problem can be addressed with differ- 
ent methods. One possibility, a representative of which 
is the Dyer-Roeder approximation [28', is to construct a 
general distance-redshift relation which would take into 
account the effect of inhomogeneities in some average 
way. However, such approaches are in general difficult to 
control |18| because they rely on approximations whose 
domain of applicability is unknown. An alternative pos- 
sibility consists in constructing inhomogeneous cosmo- 
logical models, with a discrete distribution of matter, 
and studying the impact on light propagation. Several 
models exist in the litterature: the Schwarzschild-cell 
method [IS] or the lattice universe [3D] which are both 
approximate solutions of the Einstein's equations; and 
the Swiss-cheese models [SI] which are constructed by 
matching together patches of exact solutions of the Ein- 
stein's equations. This last approach is the one that we 
shall follow in the present letter. 

Consider a Swiss-cheese model in which clumps of mat- 
ter, each of them lying at the center of a spherical void. 



are embedded in a FL spacetime. The interior region 
of a void is thus described by the Kottler geometry (i.e. 
Schwarzschild with a cosmological constant), while the 
exterior geometry is the FL one. By construction, such 
inhomogeneities do not modify the expansion dynamics 
of the embedding FL universe, thus avoiding any discus- 
sion regarding backreaction. The resuhing spacetime is 
well-defined, because the Darmois-Israel junction condi- 
tions are satisfied on the boundary of every void. Com- 
pared to a striclty homogeneous universe, a Swiss-cheese 
model is therefore characterized by two additional pa- 
rameters: the size of the voids (or equivalently the mass 
M of their central bodies), and the volumic fraction of 
the remaining FL regions which encodes the smoothness 
of the distribution of matter. It is naturally quantified 

by 
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called smoothness parameter, where Vph is the volume oc- 
cupied by the FL region within a volume V of the Swiss- 
cheese. With the definition (IH]), / = 1 corresponds to a 
model with no hole - i.e. a FL universe - while f — 
corresponds to the case where matter is exclusively under 
the form of clumps. 

The propagation of light in this model has been com- 
prehensively investigated in Ref. [23] as well as the effects 
on the Hubble diagram, generalizing earlier works |32j . 
Note that the case of voids with a Lemaitre-Tolman- 
Bondi geometry has also been considered in the litter- 
ature, but they were shown to have only a marginal ef- 
fect on the Hubble diagram [321 (and it does not question 
the fluid limit). Of course Swiss-cheese models cannot be 
considered as realistic cosmological models, but neither 
do the exact FL geometry used to interpret the Hubble 
diagram. Both spacetimes describe a spatially statisti- 
cally homogeneous and isotropic universe, and the former 
allows additionally to investigate the effect of a discrete 
distribution of matter. Since the FL universe is a partic- 
ular Swiss-cheese model, this family of spacetimes there- 
fore allows to estimate how good the hypothesis of strict 
spatial homogeneity (with a continuous matter distribu- 
tion at all scales) is. 

In our previous analysis [23], we integrated the Sachs 
equation (SJ) and simulated Hubble diagrams for a series 
of Swiss-cheese models with various cosmological param- 
eters (fimOj^Ao) characterizing the FL region, and for 
different smoothness parameters /. This allowed us to 
determine the typical error on the parameters induced by 
wrongly assuming that light propagates in a FL space- 
time, when the observations are interpreted. The effect 
of inhomogeneities is of order of a few percents, which 
can be significant given the precision on the cosmological 
parameters reached today. Such an order of magnitude is 
comparable with other estimates in similar contexts |34j . 

The synthetic Hubble diagrams enabled us to infer a 
phenomenological expression for the luminosity distance 
_Dl(z; ^mOj ^AOj Hq, f) depending on the smoothness pa- 



rameter /. This expression was then used to fit the Hub- 
ble diagram constructed from the SNLS 3 catalog [12]. 
Fig. 25 of Ref. [23] shows that / influences the result of 
the best-fit on rimo that can shift from 0.22 for / = 1 (in 
agreement with the standard FL analysis performed in 
Ref. [12]) to 0.3 for / = 0. 

We present here the constraints in the plane {Hq, fimo) 
for different values of / (assumed to be constant). The 
result is depicted on Fig. [T] which clearly shows that the 
agreement between the CMB and the Hubble diagram is 
improved for small values of /. We also note that Hq 
decreases with /, but too slightly to completely reconcile 
all the observations. 
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FIG. 1: Comparison of the constraints obtained by Planck 
on {Qn^OjHo) [Sj and from the analysis of the Hubble diagram 
constructed from the SNLS 3 catalog. The different contour 
plots correspond to different smoothness parameters /. For 
/ = 1, the geometry used to fit the data is Friedmannian, as 
assumed in standard analysis. 

In principle, the smoothness parameter / should be 
related to the growth of structures, so that we expect it 
to be an increasing function of z. A generic form 
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which interpolates between /o today and / = 1 at high 
redshift, can be used to take this effect into account. 
In this case, the impact of the inhomogeneities on 
the Hubble diagram is weaker. A redshift-dependent 
smoothness therefore reduces the agreement between 
the Hubble diagram and the CMB compared to the 
case / = /o, or alternatively requires a lower value for /q. 

Our analysis, though relying on a particular class of 
models, indicates that the FL geometry is probably too 



simplistic to describe the universe for certain types of ob- 
servations, given the accuracy reached today. In the end, 
a single metric may not be sufficient to describe all the 
cosmological observations, just as Liliputians and Broh- 
dingnag^s giants ^35] cannot use a map with the same 
resolution to travel. A better cosmological model prob- 
ably requires an atlas of maps with various smoothing 
scales, determined by the observations at hand. 

Other observations, such as lensing [Mj, may help to 
characterize the distribution and geometry of voids [37] , 
in order to construct a better geometrical model. For 
the first time, the standard FL background geometry 



may be showing its limits to interpret the cosmological 
data with the accuracy they require. 
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